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Abstract. Spin effects in a normal two-dimensional (2D) electron gas in lateral 
contact with a 2D region with spin-orbit interaction are studied. The peculiarity of this 
system is the presence of spin-dependent scattering of electrons from the interface. This 
results in an equilibrium edge spin current and nontrivial spin responses to a particle 
current. We investigate the spatial distribution of the spin currents and spin density 
under non-equilibrium conditions caused by a ballistic electron current flowing normal 
or parallel to the interface. The parallel electron current is found to generate the spin 
density near the interface and to change the edge spin current. The perpendicular 
electron current changes the edge spin current proportionally to the electron current 
and produces a bulk spin current penetrating deep into the normal region. This spin 
current has two components, one of which is directed normal to the interface and 
polarized parallel to it, and the second is parallel to the interface and is polarized in 
the plane perpendicular to the contact line. Both spin currents have a high degree of 
polarization 40 — 60%). 
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1. Introduction 

Spin transport in semiconductor nanostructures with spin-orbit interaction (SOI) 
attracts much attention, stimulated by the interesting spin dynamics in systems 
with non- conserving spin and possible applications in creating and manipulating spin 
polarization in nonmagnetic systems by purely electric methods using SOI [H Ej |3l H] . 
Spin currents are widely studied in the cases where there is a spin Hall effect, either 
intrinsic [5], El [3, El E] or extrinsic [inilllllZ]- The basic experimentally observable effect 
is the spin accumulation near lateral edges of the sample if the electric current is driven 
along it. The spin density is accumulated because the transverse spin current does not 
pass through the edges. 

A somewhat different effect appears due to electron scattering on the sample 
boundary if the particle current flows parallel to it. In this case a spin polarization 
is formed near the boundary and the spin density distribution has a rather complicated 
form [I21II31[T11[I5]. 

In this paper another system is studied, in which electrons can pass through the 
boundary of the SOI region. We consider a 2D system composed of an electron gas 
with SOI and normal (N) 2D electron gas without SOI, lying in one plane. This system 
is designated as SOI/N contact. The main purpose is to find out whether the spin 
currents, which exist in the SOI region due to the coupling between spin and orbital 
degrees of freedom, penetrate into the N region and how the spin density is distributed 
there under the non-equilibrium conditions when an electron current flows. 

This system is also interesting because it reveals a new aspect of the problem of 
equilibrium spin currents. The existence of background spin currents in thermodynamic 
equilibrium was pointed out by Rashba [16] for infinite 2D electron gas. Usaj et al [T2] 
showed that equilibrium spin current exists also in a bounded 2D electron system with 
SOI, the spin current flowing near the edges of the sample. 

The phenomenon of equilibrium spin currents has provoked an intense discussion 
on the discerning of spin currents and their definition. The question arose since spin 
currents in the SOI medium are not unambiguously defined [HI [TTl [18]. Sun et al [18] 
found that the equilibrium spin current exists in a quantum ring even if the SOI is 
absent in some segment. In this connection, of more interest is a system, in which 2D 
SOI medium is in contact with a bulky normal 2D electron gas where the spin current is 
well defined. The question is whether the equilibrium spin current exceeds the boundary 
of the SOI medium, entering into the N electron gas and what spin effects appear there. 

This question was raised in our previous papers [191 EH] . We have established that 
the equilibrium spin current does not pass through the boundary, but an equilibrium 
edge spin current flows near the boundary in both SOI and N regions. Its density 
considerably exceeds the density of the bulk equilibrium spin current. An important 
point is that the spin current is well defined in the N electron gas and its definition is 
beyond doubt. In the present work we study spin effects resulting from the interface 
scattering under non-equilibrium conditions appearing when an electron current flows in 
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Figure 1. Energy diagram of the SOI/N contact. Uq is the potential step at the 
interface; is the equihbrium chemical potential; A/i^. = — /i_ is the difference 
of chemical potentials of the right- and left- moving electrons in the presence of a 
perpendicular electron current. 



the system. This question is important since only non-equihbrium transport produces 
a spin accumulation that can be detected experimentally [Tf |2T]. 

The calculations are carried out for a weak deviation from the equilibrium caused 
by a ballistic electron current flowing normal or parallel to the interface. The electron 
current is supposed to be created by a difference of chemical potentials A/i of electrons 
moving in opposite directions. Two main effects are found. The electron current parallel 
to interface leads to the spin accumulation near the interface and changes in the edge spin 
current. The accumulated spin density is proportional to A/i, while the augmentation of 
the edge spin current is proportional to AyU^. This means that the non-equilibrium spin 
current does not depend on the polarity of the electron current. If the electron current 
is perpendicular to the interface, the edge spin currents vary linearly with A/i. However, 
the main effect consists in an appearance of the bulk spin current in N region. It has two 
components. One component is directed parallel to the interface and polarized in the 
plane perpendicular to it. The second component is directed normal to the boundary 
and polarized parallel to it. As a matter of fact, it is a spin Hall current caused by the 
SOI/N interface scattering. Both these spin currents penetrate deep into the normal 
region and have considerable polarization. 



2. The model 



Consider a 2D SOI/N contact. The energy diagram is shown in Fig. [H The SOI 
and N regions are located correspondingly at x < and x > 0. To be specific we 
suppose that the potential of the SOI region is higher than that of the N region. This 
arrangement allows one to consider an interesting case where free electrons are present 
only in the N region, though they penetrate into the SOI barrier. We assume also 
that the external electric field is too weak to change the equilibrium electron states and 
the electron current is caused by a non-equilibrium occupation of unperturbed electron 
states. Such a situation can be realized, for example, if electrons move ballistically and 
their distribution function is formed by adiabatic contacts, with the applied voltage 
dropping across the near-contact regions. 
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The wave functions are determined by the Hamiltonian: 

p"^ + p"^ a 

H = -^^^ + -{Pyax-pxCTy) + U{x), (1) 

where px,y are electron momentum components; m is the effective mass, which is 
supposed to be constant all over the structure; a^^y are the Pauli matrices; a is the 
parameter of SOI (a = for x > and a = const for x < 0); f/(x) is the potential 
energy, which is also a step function: U{x) = —UqQ^x). 

The wave functions of electrons moving to the right (R) and to the left (L) are 



{I A, kx, ky) + ^ '"a^a'I'^ ' ^a;; ky) , X < 
A' 

{\K-qx,ky) + T.rx,y\^'^1x,ky), X > 
Etj:y\x',-kxX), ^<o. 
A' 

Here kx^y and Qx are components of the electron wave vectors in the SOI and N regions, 
the index A numbers the branches of the energy spectrum and spin states. 

The reflection and transmission matrices, r^f, and t^'^,, are found using the 
following boundary conditions [22j: 



hO 



0, ^ 
dx 



^ - {iksoCTy - Pkya^)ip 



(4) 



J -0 

where fc^o = ma/ff' is the characteristic wave vector of the SOI. The parameter f3 
is introduced to take into account the SOI resulting from a potential gradient at 
the interface. This parameter is usually small. Its value is estimated as [221 120] 
f3 = 2ksoU/{eFz), with being the electric field normal to the 2D layer. For a 
typical experimental situation of InGaAs quantum wells [23] with kgo = 10^ cm~^, 
U = 5- 10-2 eV and F, = 10^ Vcm'^ one finds p = 0.01. 

The basic states \X,kx,ky) are characterized by the energy E, the wave vector 
component ky, which is invariant since the structure is uniform along the y axis, and 
the spin index A. These states were in detail described previously [191122]. Therefore, we 
only mention that there are three types of states. For E < —Ego, where Ego = ma^ / {2h^) 
is the characteristic energy of SOI, the electron states are evanescent. They are 
characterized by the complex wave vector fc^. For E > —Ego all states are either 
propagating, with real kx, or entirely decaying, with imaginary kx- Each of these states 
exists in a corresponding interval of the tangential momentum ky. If E' > 0, the index 
A = ± stands for the chirality. In an interval —Ego < E < 0, the electron states have 
only negative chirality, but there are two kind of states with a given energy. They differ 
in the wave vector kx and for convenience are also labeled by the same index A = ±. Of 
course, all these states are used in calculating the spin density and the currents. 
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Using the wave functions ([2]) and ([3]), we find the spin current density 27 and the 
spin density 5"^ for a state Iv^^'^g)? with r = {R,L) denoting the right- or left- moving 



states: 

h 



IJ{r, A, k; x) = '-^ { ^^"l a^Vi + Via^ ip^'i ) , (5) 



4 \ A,fe 



(r) 



A,fc 



h 



S^{r,X,k;x) = 'iU:la,^':i), (6) 



2 V X,k ' ^X,k/ 

where 7 labels the spin components, v is the velocity. 

The total spin current and spin density S are found by summing equations ([5]) 
and IQ over all states occupied according to the distribution function. 

In the presence of a charge current J, electrons have an asymmetric momentum 
distribution function f\{k). We simulate this asymmetry by two semicircles in fc-space 
with different Fermi momenta and, accordingly, with different chemical potentials /i+ 
and /x_ for electrons moving in opposite directions. This distribution function can 
be realized when electrons move ballistically between two contacts with the potential 
difference A/i = /i+ — In what follows we consider two cases where the electron 
current is normal or parallel to the boundary. Accordingly, all quantities (A/ij, Sj and 
J'j) are labeled by index j = {x,y) specifying the direction of the electron current. 

As a result we obtain the following expressions for the spin current and spin density: 

A/i„ x)=Y,Y. dE dkyX 

X f,{k,Afi^)D^;:\E,ky)I^{r,X,k), (7) 



poo r 
r A -^-^0 J 



X 



X h{k,Afi,)D^{\E,ky)S'^{r,X,k), (8) 

with D^[\e, ky) the density of states with chirality A. 

For small A/ij <^ Ego, Uq the spin current and density can be expanded in A^j-. 

+ gj/i-5,.,)Af^, + glfAf,l (9) 

S]ix) = x]il - - 5,A.)A/i,. (10) 

Here the expansion of J7 includes the equilibrium spin current J^^^^^ ■ A quadratic term 
is also important because some linear terms are equal to zero. To be exact, the linear 
response of the tangential spin current to the parallel current of particles is absent. 

In section 3 we study in detail the response functions of the spin density, xj; 
the spin current, Q^j. They are functions of the coordinate x. The calculations are 
carried out using the following distribution function: 



Spin currents in a normal two-dimensional electron gas in contact 



6 



Thus, the calculation scheme is the following. First, the matrices of reflection 
and transmission coefficients, r^'y and t^'y, are found from equations ([2])-(jl]). Then 
the partial spin current and spin density components, 27 ^"'^ determined using 
equations and IQ. Finally, the total spin current and spin density, J^i'jix) and Sj{x), 
as functions of x are obtained from equations (JTj) and (IHl) using numerical integration. 

3. Results and discussion 

Below, the results of calculations of the spin density and spin currents are presented for 
two directions of the electron current: normal and parallel to the SOI/N interface. 

But we begin with recalling the main results concerning the equilibrium spin 
current. 

In unrestricted 2D electron gas with SOI the equilibrium spin current has two 
nonzero components: J'y^^^ = —J^x^^^ = (2/37r)j77v, where J'n = fii^klg/Am [T6] . 
Near the SOI/N interface the perpendicular component of the spin current vanishes 
and it is absent in the N region: J^^^\x > 0) = 0. However, near to the interface 
there is an edge equilibrium spin current, which flows parallel to the boundary in both 
regions [I9]. The density of this current is considerably higher than the bulk density 
(2/37r)j77v, since the edge current is proportional to the SOI constant a, while the bulk 
spin current is proportional to a^. The spin of the edge spin current is polarized in a 
plane perpendicular to the interface and turns in this plane with increasing distance x. 
Therefore, only two components of the spin current, Jy^^\x) and Jy^^\x), are presented 
in Fig El 

The signs of the spin currents Jy{x) in SOI and N regions are seen to be opposite, 
but the total current is nonzero. It is interesting that the equilibrium edge spin current 
exists at < —Ego, i-e. even if free electrons are absent in the bulk of the SOI region. 
The current is caused by electrons penetrating under the barrier into the SOI region. 
On increasing the Fermi level above the conduction band bottom in the SOI region, 
/i > —Ego, the spin current appears in the bulk of the SOI region. It increases with fi 
to reach a limiting value {2/?)'k)Jm at = and then does not change. In contrast, the 
edge spin current increases monotonically with /i, approximately as {jj, + UoY^"^. 

Note that the parameter j3 does not affect the results very greatly. Thus, with 
varying (3 the spin current changes by an amount of the order of j3. 

3.1. Parallel electron current 

If the electron current flows parallel to the SOI/N interface, the picture of the 
distribution of edge spin currents changes quite a bit. The point is that in this case the 
linear response of the spin current is absent, Qy = Gy = 0. This occurs for the following 
reason. The edge spin current is created equally by electrons moving up and down the 
y axis. If the distribution function f{ky) is symmetric, the charge currents of the up- 
and down-moving electrons cancel each other, while the spin currents, in contrast, are 
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summed. Therefore, in the thermodynamic equihbrium the particle current is absent 
but the edge spin current is nonzero. A small asymmetric variation of the distribution 
function Afiky) results in the parallel particle current (linear in Af{ky)), but in the 
first approximation does not change the spin current. The variation of the edge spin 
current is quadratic in Af{ky). Hence it is proportional to the square of the electron 
current and in this approximation does not depend on its polarity. 

The augmentation of the spin current, caused by the parallel electron current, is 
described by the factors Qy^y'^ and Qy^y^ in equation (Q. The calculation shows that 
these factors decrease with distance from the interface and oscillate with the period 
~ Tr/kso- 

Another stronger effect consists in the spin accumulation near the boundary in both 
SOI and N regions. The spin density has two nonzero components, Sy and Sy, which 
are proportional to A^y in the first approximation. The spatial distribution of the spin 
susceptibility components Xy and is presented in Fig. [3] for different positions of the 
Fermi level. 

The spin susceptibility x^i^) is seen to oscillate and fall to zero with the distance 
X in the N region for all Fermi levels. In the SOI region, Xy{^) ^ilso drops to zero 
at —Uq < h < —Eso- However, at /i > —Ego the spin density in the bulk reaches 
asymptotically a constant value, which is determined by the difference in occupation 
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between electron states with velocities directed along y axis and opposite to it. 

The spin susceptibility component Xy{^) oscillates with a period ~ n/kso and slowly 
goes to zero with the distance in the bulk of both SOI and N regions. As a function 
of the Fermi energy at a given Afiy, the spin susceptibilities Xy and Xy increase in the 
range < —Ego and decrease at > —Ego- 

Numerical estimation of the spin density accumulated near the boundary in InGaAs 
quantum well for Uq = 6Eso, = Eso and the electric current density 1 mA/cm give a 
value of the order of 2iS/^ ~ l/im~^. Recent experiments have shown that such a value 
of spin density can be detected by Kerr rotation microscopy [3 El] . 

3.2. Normal electron current 

The electron current can flow perpendicularly through the interface only if /i > —Ego, 
i.e. when free electrons are present in the bulk of the SOI region. First, note that in 
the bulk of the SOI region there is an equilibrium spin current Jx'^^^ |I6] flowing normal 
to the interface, but this current does not pass into the N region [19]. In the presence 
of a non-equilibrium electron current, this spin current keeps unchanged if /x > 0, and 
slightly diminishes proportionally to [eV^/(/i + EsoYf' -C 1 otherwise. Thus spin effects 
in the N region originate solely from the spin-dependent scattering of electrons from the 
interface. The main effect is the appearance of the spin current in the N region. There 




Figure 4. Spatial dependence of the normal components of the spin current J^^. for 
different polarities of the electron current. Curve 1 is the equilibrium spin current 
(A/ij; = 0); curve 2 is the spin current in the presence of the electron current flowing 
from N to SOI region (A/x^^ = 2Eso)', curve 3 is spin current for the electron current of 
opposed direction {Afi^ = —2Eso)- Dashed line indicates the asymptotic value of the 
spin current in the bulk, /i = 2Eso, Uq — QEso and /3 = 0.01. 



are two components of the non-equilibrium spin current: normal and tangential. 

The spin current component directed normally and polarized in y direction, J^^, 
arises in both SOI and N regions as a result of the interface scattering. This is a 
longitudinal spin current effect. Its dependence on the distance x is shown in Fig. H] 
for two opposed directions of the electron current (A/z > and A/z < 0). In the SOI 
region the spatial dependence of J^^^ reflects the interference pattern of the incident and 
reflected electron flows. In contrast, in the N region the spin current does not depend on 
the distance, since the spin is a conserved quantity there. The spin density is constant 
as well. The spin polarization (defined as IJl^xl i^'J)) depends on the chemical potential 
reaching a maximum (about 50%) when /x lies close to the conduction band bottom in 
the SOI region. 

The tangential spin-current has two components of spin polarization lying in the 
plane perpendicular to the interface, J7^^ and J'^ .^. This is a transverse spin current 
effect. The dependence of the transversal currents on the distance is shown in Fig. [5] 
for opposed polarities of the electron current. Changes of the transverse spin current in 
the N region under the action of an electron current Jx are, of course, a manifestation 
of the spin Hall effect caused by the spin-dependent scattering from the interface. 

The mechanism of the non-equilibrium spin effects in the bulk of the N region can 
be understood in a way similar to that suggested in Ref. [19] to explain the equilibrium 
edge spin currents. Consider an unpolarized electron flow falling on the interface from 
the N region. During the scattering the electrons partially penetrate into the SOI region 
where an effective magnetic field forces their spin to precess. In the case of the Rashba 
SOI, this field is perpendicular to the wave vector and lies in the plane. As a result of the 
spin precession, the spin polarization of reflected electrons changes. Spin components 
Sx and Sz of electrons with equal kx but opposite ky turn out to be of opposite sign 




Figure 5. Distance distribution of the tangential spin currents with spin polarized in 
X and z directions (panels (a) and (b)), for different polarities of the electron current. 
Curves 1 depict the equilibrium current at A/i = 0; curves 2 and 3 are the spin currents 
for opposed directions of the electron current, = ±2Eso- Uq — 6Eso and /3 = 0.01. 



and cancel each other. Therefore under the equihbrium conditions the spin density is 
absent, but the edge spin current exists since the partial spin currents of electrons with 
opposed ky are summed up. However, in the bulk of the N region, the spin current is 
absent because the spin currents of electrons with opposed cancel each other. This 
compensation is disturbed when an electron current flows through the interface resulting 
in the net total spin current in the bulk. 

The fact that this current exists in the N region far away from the interface, where 
no SOI is present, is caused by the infinite size of the interface scatterer in the y direction. 
It is obvious, that the transverse spin currents should decay with distance deep into the 
N region because of scattering processes not considered here. 

The spin currents J'y^ and J'y^ increase linearly with A/x^. and therefore can be 
characterized by spin Hall conductivities Q^^ and Q^^. As distinct from the data of 
Fig. E], they describe only the non-equilibrium part of the spin currents. The spatial 
distributions of the spin Hall conductivities are shown in Fig. [61 

It is worth emphasizing that the spin-Hall current in the bulk of the N region 
increases with the Fermi energy even if ^ Ego approximately as /i^''^, in contrast to 
the case for the normal spin current J']^^ and the spin density S^. Numerically for the 
data of Fig. [6] the density of the transverse spin current is ~ 0.45(^/2e) J at /x = 0.5Eso 
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Figure 6. Spatial distributions of x and z components of Hall spin conductivity (a) 
and (b) for different Fermi levels /x in the presence of a perpendicular current. Curve 
1 - = O.BEso, curve 2 - /x = 2.5Eso- Uq = GE^o and /3 = 0.01. 



and exceeds 0.6(/i/2e)J at /i = 2.5-E'so- Thus this effect is more robust and we beheve 
that it could be attractive for experimental observation. The mechanism of increasing 
the transverse spin current with /i is explained taking into account that this current 
is created mainly by electrons moving nearly parallel to the interface with the velocity 
~ hkp/m. 

The transverse spin current in the N region could be investigated experimentally 
by measuring the spin density accumulated near the lateral edges of the sample. The 
density of spins near edge is estimated as follows 

with Dg the spin diffusion coefficient, Lsr the spin relaxation length. For Qy^ = O.Obkgo, 
kso = 10^ cm"\ Lsr = 2 ■ 10~^ cm, Dg = 10^ cm^s~^ and Afi = 10~^ eV one finds 
Ug = 3 ■ lO^cm"^ = 30/im~^. This value of spin density is readily measurable in the 
experiments (TJ [21] . 



4. Conclusion 



We have found substantial spin effects in a normal 2D electron gas in contact with a 
2D medium with SOI, which lies in the same plane. The spin effects are caused by the 
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spin-dependent scattering of electrons from the interface. 

In the thermodynamic equihbrium state, the boundary scattering results in the 
appearance of an edge spin current flowing along the boundary in both SOI and 
N regions, with spin polarization being directed in the plane perpendicular to the 
boundary. 

Under the non-equilibrium conditions caused by a ballistic electron current flowing 
in the system perpendicularly or tangentially to the interface, the following effects are 
developed. The tangential electron current gives rise to the formation of the edge spin 
density in both SOI and N regions. In addition it leads to the non-equilibrium edge spin 
currents that are proportional to the square of the electron current. 

If the electron current flows normal to the interface, the main effect is the 
appearance of the non-equilibrium spin current. It is proportional to the charge current. 
The spin current has two components: longitudinal and transverse with respect to the 
charge current. The longitudinal component is polarized parallel to the boundary. The 
transverse component is polarized in the plane perpendicular to the boundary. By its 
origin this component is the spin Hall current caused by the spin-dependent scattering of 
electrons from the interface. In view of possible experimental realizations, it is essential 
that the spin currents penetrate into the N region to a distance which is determined 
by scattering processes in the bulk and can really reach several microns. The spin- 
polarization degree of the currents is high enough (about 40 — 60%). 
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